ON THE SELF-DECOMPOS ABILITY OF THE FRECHET 

DISTRIBUTION 



PIERRE BOSCH AND THOMAS SIMON 

Abstract. Let {T t , t > 0} be the Gamma subordinator. Using a moment identification 
due to Bertoin-Yor (2002), we observe that for every t > and a € (0, 1) the random variable 
IY" is distributed as the exponential functional of some spectrally negative Levy process. 
This entails that all size-biased samplings of Frechet distributions are self-decomposable and 
that the extreme value distribution is infinitely divisible if and only if £ ^ (0, 1), solving 
problems raised by Steutel (1973) and Bondesson (1992). We also review different analytical 
and probabilistic interpretations of the infinite divisibility of r^~ a for t, a > 0. 



1. Introduction 

The extreme value theorem - see e.g. Theorem 8.13.1 in jl] - states that non-degenerate 
distribution functions arising as limits of properly renormalized running maxima of i.i.d. 
random variables belong to one of the families 

F (x) = e~ e ~\ xGl, or F e (x) = [ 1 ,% ^ ^ > ° , x > 0. 
uv ' ' ' I e-' x 75 if £ < 

The distribution Fq is known as the Gumbel distribution, whereas F^ is called a Weibull 
distribution for £ > and a Frechet distribution for £ < 0. In the following, we denote by 
X% the random variable with distribution function F%. Observe that 

^- ~ d v „ 
— 7— 1 — > X Q as £ -> 0, 

so that the above parametrization is continuous in £. In the present paper we are interested 
in the self-decomposability (SD) of X^, referring e.g. to Section 15 in [2] for an account on 
self-decomposability. The Gumbel distribution is SD because of the identities 

X Q = — log L = —a log L + a log S a 

for every a £ (0,1), where here and throughout L stands for the standard exponential 
variable and S a for the standard positive a-stable variable - see e.g. Exercise 29.16 in [T3] 
for a proof of the second identity. If £ £ (0, 1) then the variable X^ is not infinitely divisible 
(ID) and hence not SD, because of its superexponential distribution tails - see e.g. Theorem 
26.1 in P3]. When £ > 1, the variable X^ has a completely monotone density and is ID by 
Goldie's criterion - see e.g. Theorem 4.2 in [17], or by the ME property which makes it the 
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first-passage time of some continuous time Markov chain - see e.g. Chapter 9 in [5] for an 
account. When |£| > 1, the identity in law 

X 6 = L« 

and the HCM theory of Thorin and Bondesson [5j show that the distribution of is a 
generalized Gamma convolution (GGC) and is hence SD - see Example 4.3.4 in [5]. The 
natural question whether X% is SD or even ID for £ G (—1,0) was first raised by Steutel in 
1973 - see Section 3.4 in [T7], and has remained open ever since. In section 4.5 of [5] - see 
also the Appendix B.3 of [18], this problem is rephrased in the broader context of generalized 
Gamma distributions. The latter are power transformations of T t where {T t , t > 0} is the 
Gamma subordinator, and can be thought of as size-biased samplings of X^ when £ < 0, in 
view of the formulae 

, E[f(XAX?] 

nmn - 

valid for every / bounded continuous and t > 0, with u = (t — l)/£. Recall in passing that 
Steutel's equation - see e.g. Theorem 51.1 in [T3j - establishes a precise link between size- 
biased sampling of order one and infinite divisibility for integrable positive random variables. 
In this note, we provide an answer to the above questions of [TFJ [5]. 

Theorem. For every £ G (—1,0) and t > 0, the random variable rf is SD. 

As a direct consequence of this result, all Frechet distributions are SD and the extreme 
value distribution is ID if and only if £ ^ (0, 1). Contrary to the case |£| > 1, our argument 
is probabilistic and consists in showing that rf is distributed as the exponential functional 
of some spectrally negative Levy process. This extends a classical result of Dufresne [6] for 
the case £ = — 1. The identification is made possible thanks to a entire moment method due 
to Bertoin-Yor [3], which applies in our context study. The proof is given in the 

next section. 

In Section 3, we review the possible interpretations of the infinite divisibility of rf for £ < 0. 
The classical case £ = — 1 allows at least four different formulations in terms of processes, and 
also an explicit computation of the Levy density which shows the GGC property without the 
HCM argument. For £ < — 1 the ID property is only known by analytical means and there is 
no direct probabilistic explanation, save for the case t — 1 by subordination or, tentatively, 
the spectral theory of a certain spectrally positive Markov processes. The situation for 
£ G (—1, 0) is exactly the opposite since in addition to the exponential functional argument, 
the ID property can also be obtained rigorously by a first-passage time argument for a 
spectrally positive Markov processes. On the other hand there is no analytic proof of the ID 
property for £ G (—1, 0). In this situation the GGC character of rf remains in particular an 
open question, which we plan to tackle in some further research. 

2. Proof of the Theorem 
We begin with a computation on the Gamma function. 
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Lemma. For every a 6 (0, 1) and u, t > one has 

where 

e^ a >(a + e x ' a + t(l - e x ' a )) 
ta ^ X) ~ ar(l-a)(l-e-/°)-+ 2 
is £/ie density of a Levy measure on (— oo, 0). 

Proof. We set A = t + au > and compute 

T(A + a) \p(\ + a,l-a) 



T(A) r(l-a) 

A /'+ 00 e -{a+X)x 



dx 



T(l-a)J {l-e~ x r 

'l-e-^ —dx 



a 



T(l-a)J v J (l-e~ x ) 

where the second equality comes from a change of variable and the third from an integration 
by parts. This yields 



uT{t + a(u + l)) au f +co „-tt+ au )x^ e 



-cm: 



1 _ e ~{t+w)x\ f . 

T{t + au) T{l-a)J x ' (1 - e -*)«+i 

r(f + a)\ an I* 1 c («+*>- 



,m + ^7 r/ (l-e aux )- r^rdx 

T(t) J T(l-a)J_J >{l-e*Y+i 

u + I (e ux - 1 - ux)f a , t {x)dx 



r(t) 

where again, the second equality comes from a change of variable and the third from an 
integration by parts. 

□ 

Remarks, (a) The above proof follows [2J p. 102. Notice in passing that some computations 
performed in [2] are slightly erroneous. For example the subordinator whose exponential 
functional is distributed as r~ a (with the notation of [2]) has no drift, but it is also killed at 
rate 1/T(1 — a). 

(b) The above decomposition extends to a = 1 since 

uF(t+(u + l)) , . 
T(t + u) 

is the Levy-Khintchine exponent of a drifted Brownian motion (the latter was already noticed 
in [5] - see Example 3 therein - in order to recover Dufresne's identity). However, such a 
formula does not seem to exist for a > 1. 
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End of the proof. Fix £ G (—1, 0), t > 0, and set a = — £ G (0, 1) for simplicity. The entire 
moments of T° are given for every n > 1 by 

T(t + an) 



E[r 



an] 



T(t) 

Ttt + at) TCt + an) . . . ip(n - 1) 

X • ■ • X ; TT = W 



r(t) r(t + a(n-l)) (n-1)! 

with the notation 

by the Lemma, and 

m = r = V (0+)- 

It is clear that ip is the Levy-Khintchine exponent of a spectrally negative Levy process Z 
with infinite variation and mean m > 0. By Proposition 2 in [3], this entails 

E[I7 n ] = E[/-"] 

for every n > 1, where / is the exponential functional of Z: 



oo 

„-z 3 



ds. 



Since Z has no positive jumps, Proposition 2 in [3] shows also that the random variable 1/1 
is moment-determinate, whence 



rt ^ /. 

The self-decomposability of / is a direct consequence of the Markov property. More precisely, 
introducing the stopping-time T y = inf{s > 0, Z s = y} for every y > 0, the fact that 
Z s — > +oo a.s. as s — > +oo and the absence of positive jumps entail that T y < +oo a.s. 
Decomposing, we get 

e~ Zs ds + / e~ Zs ds = / e~ Zs ds + e~ y / e~ zi ds 



where Z' is an independent copy of Z and the second equality follows from the Markov 
property at T y . This shows that for every c G (0, 1) there is an independent factorization 

I = cl + I c 

for some random variable J c , in other words that I — rf is self-decomposable. 

□ 

Remarks, (a) By the above Remark 1 (b), it does not seem that rf is distributed as the 
exponential functional of a Levy process for £ < —1. It would be interesting to have an 
explanation of the infinite divisibility of rf in terms of processes when £ < — 1. See next 
section for a more precise conjecture in the case t — 1. 
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(b) The self-decomposability of S% for every a G (0, 1) has been shown by Patie [12] in using 
the same kind of argument. Specifically, one can write 

E[g;°] = r(l + ») ?(!)■■ ■»(,,- 1) 



r(l + cm) (n-1)! 

where we use the same notation as above and, correcting small mistakes made in Paragraph 
3.2 of H21, 



u 







(1 - a)e x / a ((2 - a)e x ' a + (1 - e x / a )) 



= wr-, \ + / (e ux - 1 - ux) — — -— T —z dx 

yj r(l + a) 7_ 00 a 2 r(l + a)(l - e*/ a ) 3 - a 

is the Levy-Khintchine exponent of some spectrally negative Levy process with positive 
mean. Setting a = t = 1/2 and comparing the above expression to the one in the Lemma, 
one can check the well-known identity in law 

(2-1) = ~^=- 

Z V L 1/2 

The present paper shows that all positive powers of Su2 are SD and one may wonder if the 
same is true for S a with any a G (0, 1). See [9] for related results and also for a characteri- 
zation of the SD property of negative powers of S a when a < 1/2. 



3. Further remarks and open questions 

In this section we would like to review several existing or tentative approaches for the ID, 
SD and GGC properties of the distribution of rf or when £ < 0. 



3.1. The case £ = 0. This is a rather specific situation but we include it here for complete- 
ness. As mentioned in Section 3.4 of [17], the SD property of the two-sided X is a direct 
consequence of the extreme value theory because 



(3.1) 



L 1 + f + 



■ + 



L„ 



n 



logn = max(L 1; . . . , L ri ) — logn 



as n — > +oo, where Li,...,L n are independent copies of L ~ Exp (1). The above identity 
and convergence in law, readily obtained in comparing Laplace transforms and distribu- 
tion functions, yield after some further computations the following closed expression for the 
Laplace transform of Xq\ 



T(l + A) = exp 



- 7 A + / ( 



-Xx 



1 + Xx)- 



dx 



x\& 



"1)J 



where 7 is Euler's constant. The complete monotonicity of l/(e x — 1) shows then that Xo 
is an extended GGC in the sense of Chapter 7 in [5]. See also Exercise 18.19 in [14] and 
Example 7.2.3 in [5] for another argument based on Pick functions, recovering (13. ip . 
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3.2. The case £ = — 1. This is the classical situation, very well-known, but we give some 
details for comparison purposes. The ID property of can first be understood by the sole 
fact that 

i im (j?E—\ 



n->+oo y 1 + nX 

because the left-hand side is the first-passage time distribution function of a certain birth 
and death process - see Theorem 3.1 and (3.3) in [17]. The random variable T^ 1 is also a 
GIG and is hence distributed as the unilateral first-passage time of a diffusion [T], which 
explains its infinite divisibility by continuity and the Markov property. More precisely one 
has 

(3.2) ±- £ mi{u > 0, XI = 0} 

where {X l ul u > 0} is a Bessel process of dimension 2(1 — t) starting from one. The SD 
property follows as above from Dufresne's identity [6], which reads 



CO 



(3.3) 1 = / e B -- tu ' 2 du 

M Jo 

where {B u , u > 0} is a standard linear Brownian motion. Also, Exercise 16.4 in [H] shows 
that T^ 1 is the one-dimensional marginal of a certain self-similar additive process, whence 
its self-decomposability by Theorem 16.1 in [H]. The link between this latter interpretation 
and f)3.2p and (13. 3 j) has been explained in depth in [19]. 

It does not seem that these four interpretations can provide any explicit information on 
the Levy-Khintchine exponent of I^T 1 . But in this case analytical computations can also 
be performed. More precisely, taking for simplicity the same normalization as in (I3.2p and 
setting ipt(X) = -logE[e- A / 4rt ], formulae (7.12.23), (7.11.25) and (7.11.26) in [7J entail 

(3.4) m- KM 



2VXK t (VX) 



where K t is the Macdonald function. This shows ^/ 2 (A) = 1/2VA viz. v 9 i/2(A) = vA when 
t = 1/2, and one recovers the identity (12.11) . For t = 3/2, one obtains 

^ 3/2K ' 1 + L J J \X + xJ n(l + x) 

where the first equality follows from Formula (7.2.40) in [7J, and the third equality from 
Exercise 29.16 in [H] and (2.2.5) in [5]. This means precisely - see (3.1.1) in [5] - that the 
distribution of is a GGC with zero drift and Thorin measure 



, . dx 
U ^ {dx) = 2rtT+xj' 

The latter property can be extended to all values of t thanks to a result originally due to 
Grosswald [8] on Student's distribution. Together with f l3.4jl . the main theorem in [8] entails 
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namely that the distribution of l/4T t is a GGC with zero drift and Thorin measure 

1 



U t (dx) 



where Jt and Y t are the usual Bessel functions of the first kind - see [7] p. 4. 

3.3. The case £ G (—1,0). In this situation, the present paper yields an interpretation of 
the self-decomposability of rf by the identification 

POO 

rf = J e- z «du, 

where Z is a spectrally negative Levy process. Another explanation, similar to (13. 2p . can 
then be obtained by the Lamperti transformation - see e.g. the introduction of j3] for an 
account and references. More precisely, setting 

Y u = exp[-Z Tu ] 

with the notation t u = inf{s > 0, J Q S e~ Zv dv > u}, then Y = [Y u , < u < rf} is a 
spectrally positive Markov process (which is also self-similar) starting from one and we have 

rf = mi{u > 0, Y u = 0}, 

so that the infinite divisibility of rf (but not, or at least not directly, its self-decomposability) 
is a consequence of the Markov property and the absence of negative jumps for Y. It would 
be interesting to see if rf could be embedded in some self-similar additive process analogous 
to the Brownian escape process of the case £ = 1, described in Exercise 16.4 of |14j . 

Our main result can also be interpreted analytically in terms of generalized Bessel func- 
tions. Setting a = — £ and writing down 

(3.5) E[e~ xr '} = -4rr / x-^e-^^dx - lM ! 



aT(t) J aT(t) 
with the notation (1.7.42) of [TT], the infinite divisibility of rf entails that the function 

'<(A)' 



(3.6) A ^ 



Z pW 



is completely monotone for any p > 1 and v > 0. One might ask if the latter function is also 
a Stieltjes transform, which is equivalent to the GGC property for the distribution of rf - 
see Chapter 3 in [5]. Indeed, it is very natural to conjecture such a property for £ G (—1, 0) 
in view of the above cases £ = and £ = — 1. Compared to classical Bessel functions, the 
theory of generalized Bessel functions is however rather incomplete, and proving like in [8] 
that the function (13. 6p is a Stieltjes transform is believed to be challenging. 
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3.4. The case £ < —1. In this situation the GGC property of the distribution of rf is 
most quickly obtained from the HCM character of the density function - see Chapter 5 and 
especially Example 5.5.2 in [5J. Notice that this analytical argument extends to £ = — 1 but 
not to £ G (—1,0) since otherwise would also have a HCM density and hence be ID, 
which is false. This entails that the function in (13.61) is indeed a Stieltjes transform for any 
p G (0, 1) and v > 0, and it would be interesting to identify the underlying Thorin measure 
as in Grosswald's theorem. 

A probabilistic interpretation of the self-decomposability of T\ = can also be given by 
Bochner's subordination. Setting a = — 1/£ G (0, 1), one has indeed 

r£ A r-1 v c A qa 

J-J — J-J A O a — O j^ — a • 

where {S 1 ", u > 0} stands for the a— stable subordinator with marginal S® — S a . Since L~ a 
is SD by our result, this means that Lfc is the marginal of some subordinator which is itself 
subordinated to the a— stable one, and Proposition 4.1. in [T5] shows that is SD. Besides, 
setting if a resp. ip^ for the Levy-Khintchine exponent of L~ a resp. one deduces from 
Theorem 30.4 in [H] the following relationship 

<^(A) = if a (X a ). 

Another, tentative, probabilistic interpretation of the self-decomposability of could be 
given in terms of a certain spectrally positive Markov process. Setting a = — 1/£ and 
Va(^) — E[e~ ALC ], Theorem 4.17 p. 258 in [11] and (13. 5p above show that y a is a solution to 
the fractional differential equation 

xD°: +1 y a - ay a = 0, 

where D® +1 is a fractional Riemann-Liouville derivative - see Section 2.1 in [11] . When 
a = 1 viz. £ = — 1 the above amounts to a Bessel equation and Feller's theory applies, 
making L~ x the first-passage time of a Bessel process of index - see [10]. When a G (0, 1) 
the operator _D°; +1 is the infinitesimal generator of a spectrally positive (1 + a)— stable Levy 
process reflected at its minimum, which is a spectrally positive Markov process - see Section 
3 in [13] and the references therein. By downward continuity, one may wonder if cannot 
be viewed as the first-passage time of some scale-transformation of the latter, eventhough 
no Feller's theory is available for fractional operators whose order lies in (1, 2). 

The above probabilistic interpretations do not seem to hold for t ^ 1. On the one hand, 
Theorem 4.17 in [TT] yields then an equation with two fractional derivatives of different 
order for the Laplace transform of rf. On the other hand, keeping the notation a = — 1/£, 
Theorem 1 in [16] shows the factorization 

rf = x 5(f> 

where Sa is the size-biased sampling of S a of order —at, viz. 
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The GGC character of S& follows from that of S a and Theorem 6.2.4. in [5], which shows 
by the above case £ = — 1 that Tj is the independent product of two SD random variables. 
By Theorem 16.1 in [H], this entails that there exist two independent 1-self-similar additive 
increasing processes Y and Z such that 

r? ^ y Zi . 

Unfortunately, contrary to Bochner's subordination the independent composition of two 
additive processes is not necessarily an additive process anymore, so that the above identity 
does not provide a probabilistic proof of the self-decomposability of r|. 

Acknowledgements. The authors are grateful to L. Bondesson for the interest he took 
in this work, and several useful comments. Ce travail a beneficie d'une aide de l'Agence 
Nationale de la Recherche portant la reference ANR-09-BLAN-0084-01. 
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